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Abstract 

We are used to thinking of an operator acting once, twice, and so on. 
However, an operator acting integer times can be consistently analytic con- 
tinued to an operator acting complex times. Applications: (s,r) diagrams and 
an extension of Fractional Calculus where commutativity of fractional deriva- 
tives is preserved, generating integrals and non-standard derivations of theo- 
rems in Number Theory, non-integer power series and breaking of Leibniz and 
Chain rules, pseudo-groups and symmetry deforming models in particle physics 
and cosmology, non-local effect in analytic continued matrix representations, 
particle-physics-like scatterings of zeros of analytic continued Bernoulli polyno- 
mials ( |physics / 9705021 ) , analytic continuation of operators in QM, QFT and 
Strings. 
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1 Introduction: 

Analytic continuation of Operators 



1.1 Questions in Context 

rf(l/2+i) 



What is — , ? Is it meaningful? What about / (dx) 
dxW 2+l > J 

Given a function f(x), how do we evaluate 



(1/2+0 ? 



/7(V 2 +0 , 

jrf{x) and / /(x) (dxjW^ ? 



What do they mean? 

Are these useful? Are there examples of applications? 

If A is a generic operator, how do we compute J 4( 1 / 2 +') ? 

For the creation and annihilation operators in Quantum Mechanics and 
Field Theories, how do we calculate and interpret 

a (l/2+0 \ n ) and the 

commutator [a^ 2+i \ (at) (1/2+i) ] ? 



• What about similar generalisation to other operators, eg. Supersymmet- 
ric operators, Vertex operators, Virasoro algebra in String Theories, and 
Superconformal algebra in Superstring? 

• What are their surprising implications and consequences? 
The aim of this paper is to address these issues and questions. 



1.2 The Idea 

We know that in Complex Analysis functions can be analytic continued 
from integer points n on the real line to complex plane s, eg. from fractorial re! 
to Gamma function T(s), and from within the area of convergence to beyond via 
functional equations, eg. the functional equation of the Riemann zeta function 

C( S ) = 2 s vr s - 1 sin^r(l- S )C(l- S ) Va) (1) 

We too know that in Euclidean Quantum Gravity [||, spacetime can be 
analytic continued from 

Lorentzian spacetime — > Complex spacetime — > Euclidean spacetime 

by rotating the signature of spacetime metric from 

(-1,+1,+1,+1) - (e ie ,+l,+l,+l) - (+1,+1,+1,+1) (2) 
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Path integrals ill-defined in Lorentzian spacetime become well-behaved in Eu- 
clidean spacetime in which they can be evaluated by methods of infinite-descent 
and saddle-point approximation. This gives us a tool to explore non-perturbative, 
non-linear and topological structures like instantons and wormholes, the Ther- 
modynamics of Hawking radiation ||, the end state of black hole evaporation 
Q, and Conjectures on the boundary conditions of Complex spacetime ||. 

Now, we take analytic continuation a step further. 

Can operators be analytic continued? 

Operators and their representations permeate almost every branch of Math- 
ematics and field of Sciences. If the analytic continuation of operators can be 
consistently defined and computed, then the idea may have broad applications 
and universal implications. 

We begin the first step with the two most familiar operators of all, the 

d f 
differential operator — and integral operator / dx . 

dx J 



d r 

2 Analytic continuation of — and / dx 

dx J 

2.1 Differentiating and Integrating in non-integer s-dimensions 

Analytic continuation of differentiation and integration to non-integer dimen- 
sions is straightforward. 

Differential of an integer n-dimensional function in n-dimensions is 

Q Q Q 

■ f(xi,X 2 , - ■ ■ ,Xn) (3) 



dx\ dx2 dx r 
and the corresponding integral is 

cx„ rX2 rxi 

f(x 1 ,x 2 , • • • ,x n )d n x = / ••• / / f(x 1 ,x 2 , ■ ■ ■ ,x n ) dxidx 2 ■ ■ ■ dx n 



n-times 

If / is spherically symmetric, / = f(r), then 

d d d d n ~ l d r(n/2) 



(4) 



(5) 



dx\ dx 2 dx n Q r n-i dQn-i 2-7r n / 2 dr 71 ^ 1 
can then be analytic continued to differential of a non-integer s-dimensional 
function in s-dimensions as 

5 s - 1 d r(s/2) d s ~ l 



dr*- 1 dSls-l 2tt s / 2 dr s ~ l 
and the corresponding integral 

coo r2ir 



(6) 



d n x = / r n - l dr d8i sm6 2 d9 2 --- sin"" 2 n _i^ n _i 
Jo J J J 
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oo f o_n/2 poo 

n-lj„ I ,70 . _ / „n-l 



r" ^dr / dfi n _i = — — - / r n - L dr (7) 
J T(n/2) 7o 

as 

27r s / 2 Z" 00 

" r(s/2)i (8) 

where s can be real or complex. 

An important application example of such analytic continuation is 't Hooft 
and Veltman's Dimensional Regularization where it is used to isolate singular- 
ities in divergent integrals in Quantum Field Theory ||. 

However, this is not the only possibility. 

There is another possible analytic continuation. 



2.2 Differentiating and Integrating complex s-times in one-dimension 

Think of the differential and integral as operators. Differentiating or integrating 
n-times a one- variable function f(x) can be thought of as letting the operator 
act n-times or n-fold on the function, 

d \ n d n 



(rx \n rx 

J (dx)j : f(x) ^ J f{x) {dx) n 



X fX n fX n -l fXS fX2 

f(xi) dx\ dxi ■ ■ ■ dx n -i dx n (10) 



n-times 

Note that the limits of integration of this analytic continuation are different 
from those of integrating in n-dimensions. 

• Integration in n-dimensions is a product of integrals 

whereas 

• Integrating n-times gives a set of n nested integrals with the limits taken 
at the end of the integration. 

At this point, it is natural to generalise and combine both the differential and 
integral operators into one fundamental operator 

d s 

dx^ ' Re(s)> ° 

1 , Re(s) = (11) 



D s 



J (dx)~ s , Re(s) < 



The analytic continuation of the differential and integral operators to the 
D s operator is known as Fractional Calculus 0, a subject of active research 
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and of current interest because of its widespread pure and applied applications 
in Maths, Physics, Engineering and other Sciences B. 

In this paper, an extension of the conventional Fractional Calculus is intro- 
duced. Evaluating D£x w is simple and obvious for Re(w) > . On the other 
hand, for the case of Re(u>) < , it is not so straightforward but will turn out 
to be simple when mapped to the (s, w) diagrams to be introduced. 

The surprise is that the consequences of analytic continuation of these op- 
erators are not only highly non-trivial but useful. 

In particular, a number of known and new results in Number Theory are 
derived in non-standard ways using the idea of analytic continuation of opera- 
tors in Sections ||| and ||[ These results demonstrate the usefulness and justify 
the purpose of the idea. 

In addition, when D s acts on a standard power series, the result is a non- 
integer power series. Analysis in Section [| shows that there are interesting 
relations between non-integer power series and the usual integer power series in 
limiting forms. In fact, Fractional Calculus can be reinterpreted as differential 
and integral operators acting non-integer times. D s is observed to break Leibniz 
rule and Chain rule when s is non- integer, and thus we are unable to evaluate 
directly the action of D s on a function of functions. However, by a trick of 
series expansion, we can express D s as a nested sum of d n /dx n or J (dx) n which 
we can evaluate directly with Leibniz rule and Chain rule. 

Existing concepts in Group Theory are then extended in Section ^ using 
these results. An extension of Dirac Algebra from the Dirac equation is found 
in Section [?]. In Finite Difference, the matrix representations of d n /dx n and 
/ (dx) n are sparse. When the operator D s is casted in matrix representation as 
in Section ||, the matrix becomes dense for non-integer s, and so the local finite 
difference becomes non-local. The cause of this non-local effect can be traced 
to the appearance of non-integer power series. 

Towards the end of the paper, problems were raised and challenges were 
posed on analytic continuation of operators and algebras in Quantum Mechan- 
ics, Supersymmetry, and Quantum Field and String Theories. 

All in all, the analytic continuation of operators turns out to be quite a 
general and powerful tool to explore Number Theory, Group Theory, Algebra, 
Finite Difference and Matrix Representation. Exploring with the idea has mo- 
tivated the introduction of a few other new ideas and concepts into these fields, 
each is of very different nature from the other. The intriguing results suggest 
that the idea of analytic continuation of operators may well find interesting 
widespread applications in various other fields. 
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3 Defining the action of operator Z>' 



The operator DJ} for integer n is well-defined since it corresponds to for 

dx n 



n > 1 and J {dx) n for n < -1. 



d n 
dx 11 



ml 



{m — n)\ 



x 1 



x m {dx) 



ml 



, < m < n 
, m > n > 

, m > , n < 



> , m,n £ Z 



{m—n)\ 

If we tabulate D™x m for integers n,m , we observe a pattern emerges. 



Table 1: Tabulated results of D?x r ' 









D£x 


m 












m\n 


-3 


-2 


-1 





1 


2 


3 


2 


21/5! x 5 


2!/4! x 4 


2!/3! x 3 


x^ 


2! x 


2! 







1 


1/4! x 4 


1/3! x 3 


1/2! x 2 


X 


1 















1/3! x 3 


1/2! x 2 


X 


1 















-1 


J log x {dx) ' A 


/ log x (dx) 


log X 


x- 1 


-x~ 2 


2! x" 


-3 


-3! x" 


-4 


-2 


-J log x (da;) 


-logx 




x- 2 


-2! x~ 3 


3! x 


-4 


-4! x" 


-5 


-3 


1/2! logx 


1/2! x _i 


-1/2! x" 2 


x~ 3 


-3!/2! x" 4 


4!/2! x" 


-5 


-5!/2! x" 


-6 



3.1 Fractional Calculus and its Extension 

There are several possible ways for analytic continuing d/dx and J dx : 

1. Riemann-Liouville Fractional Calculus 

In the conventional Riemann-Liouville Fractional Calculus 0j, we start 
with the Riemann-Liouville Integral 

X rx n rx n -i rx.3 rX2 \ rx j (+\ 

/ / •••/ / f {xi) dxi dx 2 ■ ■ ■ dx n -! dx n = -— 

a Ja Ja Ja >- [Jl) J a \X tj 



n-times 

(12) 

as the fundamental defining expression. 

Fractional integral is analytic continued from this Riemann-Liouville in- 
tegral (pi) as 

D^J{x) = d ,f_ a y K*) = [f{x){dx)-° 

i r fit) 



T{-a)J a Jx J lh dt f ° r ° <Q > CT ' aG R(13) 
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Fractional derivative is in turn derived from the fractional integral ( O ) 
by ordinary differentiation as (choose m> a) 

DZf{x) = D™ (D-^-^f{x)) for a > 0, m e Z+ 

and it has the property 

Lgftfx) = f3°Dfr x) f([3x) 

However, for D£x r , the Riemann-Liouville Integral definition is well- 
defined only for the half plane a £ R , r > — 1. 

2. Fractional Calculus by Cauchy Integral 

Cauchy Integral for an analytic function f(z) in the complex plane is 



f (n) (zo) 



T(l + n) 



/(*) 



1m Jc [z - z ) 1+n 



dz 



(14) 



Generalization of n to non-integer values is however not trivial as the 
term (z — zo) 1+a may become multi- valued and the result will depend on 
the choice of branch cut and integration path. 

3. Fractional Calculus by Fourier Transform or other Integral Transform 
In Fourier Transform, 



+oo 



f(x)e ikx dx , f{x) 



1 

2^ 



f{x) e~ ikx dk (15) 



-f oo 



f{k)DZ(e~ ikx )dk , (J £ R 
(-ik) a f(k) e~ ikx dk 



+oo 



{-ikfe 



-ikx 



dk 



(16) 



This was shown by Zavada |1C] to be equivalent to the Riemann-Liouville 
Fractional Calculus and the Fractional Calculus by Cauchy Integral. 

Now, to computed numerically the fractional derivative or fractional in- 



tegral of a function f(x) or multi-variable function f(x^\x^- 



we choose to use the Discrete Fast Fourier Transform (DFFT) instead. 
Take the function values in the interval or region of interest, identify the 
boundaries so as to make it periodic, feed this into the DFFT algorithm 
to Fourier transform the function into fc-space, multiply the component 
corresponding to k with (—ik) a , and feed the result into the Inverse DFFT 
to Inverse Fourier transform back it to x-space. 
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4. Extended Fractional Calculus 



Now, we introduce the Extended Fractional Calculus in which the limit 
of the ratio of Gamma functions 



T(l+r + e) 



e^o r(l+r+e-cr) 



lim 



(17) 



is taken as the fundamental defining expression instead. 

Unlike the Riemann-Liouville Integral, the limit is well-defined for the 
entire plane a,r € R except along the line intervals r £ Z~, a £ R\Z. 
The analytic continuation to these line intervals will be derived in Section 



3.2.2, equation (pq). 

However, in the wedge-shaped region r < a , r > 0, there are actually 
two possible choices, of which one (Type I) corresponds to the Riemann- 
Liouville Fractional Calculus in that region where fractional derivatives 
do not generally commute, and the other (Type II) an extension of it 
where commutativity of the fractional derivatives is preserved. This will 
be described in Section |3.2.3 . 



Fractional Calculus by Nested Series Expansion 

Here, we introduce the method of analytic continuation of operators by 
nested series expansion with the following observation: 

D s can be formally expanded into a nested series as 

1 



D s 



(wl - [ 



wl-D 



1 + 



oo 



n=l 



(-iy 



' oo 
\ n=l 



IV. 



(-iy 



in 



= w 

1-1 

IK- 

Lfc=0 

n(- 

fc=0 



■k) 



■k) 



D 

w 

Id 

w 




(18) 



where s, w £ C, and 1 is the identity operator. 

In the nested series on r.h.s, all the operators -D's are raised to inte- 
ger powers m, and D m corresponds to ordinary m-fold differentiation if 
m < or integration if m > . The region of convergence in s and the 
rate of convergence of the series will be dependent on parameter w, and 
the function on which D acts. 

In this way, Fractional Calculus with its analytic continued operator D s 
can be reinterpreted as operator D acting non-integer or complex s-times 



on a function. Examples of applications can be found in Sections 5.2 and 
0, and in g. 



3.2 Analytic continuation from integer n to real a 

To analytic continue to D°x r for real <r, r, we introduce the (a, r) diagram in 
which the coefficient of F)°x r is mapped to the point at coordinate (er, r) of the 
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diagram. The (a, r) diagram of D°x r can then be characterized into 5 regions 
as in Figure [|: 




I — F — I — I — 



Figure 1: (cr, r) diagram of D a x x r 



upper 
lower 
log 
zero 



region bounded by 



r > a , 


r 


> o 


r < a , 


r 


< 


r > cr , 


r 


< 


r < a , 


r 


> o 



A point lying on the 
3.2.1 Upper region 

r(i+r) 



right (cr > 0) is a differentiation 
left (a < 0) is an integration 



The ratio 



r(l+r-er) 
is finite everywhere in r > a , r > 
So, 

r(i+r) 



in the upper region. 



(19) 



(20) 



D?x r 

T(l + r-CT / 

3.2.2 Lower and Log regions 

Define 

f2 r <o as the union of lower and log regions such that r < , a £ R , 
^horz as the set of the horizontal lines r = — n + , n + G Z + , 



S 



VLgrid as the set of integer grid points in £l r <o such that a € Z , r € Z 

^diag as the set of right-sloping diagonals in the lower region 
r = a — n+ < , ra+ G Z + , cr, r G R . 



1 1 - : 


'—A 




' .i' ,f i ,i 4 j. 










/ / ,-' / t 4 



H>c> 



•horz 



diag 



o grid points in log region 
a grid points in lower region 



Figure 2: Lower and Log regions. 



The ratio (^) is finite everywhere in ri r< o except in ^horz\^grjd (along the 

horizontal lines "mod" the grid points). It is zero in QdiagX^grid (along the 
diagonals in the lower region "mod" the grid points). 

However, the limit 

r ( i +r+e) 

e^o T(l+r+e-cr) v y 

evaluated at the grid points in Q gr id is convergent. 
Thus, 

D°x r = lim r ^ + I +e) s ^ in fi r <o\(n ( 22 ) 

From Table |l[ the natural analytic continuation of D°x r in the part of 
tthorz\Qgrid tying in the log region is 

lim r ( 1 + r+£ ) / logxidx)^-^ = lim r ( 1 + r+e ) D (*-r) , (23) 
e^o T(e) J & { J e^o r(e) x V ; 

It can be shown using ( f22|) and the following expression 

logs = logx + log 1 = lim / x~ 1+e dx = lim - (x e — 1) (24) 

that 

J log x(dx)P = D-p log x = T ^ +p) (logx-(V(l+p)+7)) (25) 
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where the digamma function |]T 



dz ov w/ r(z) 

the Euler constant 7 = — r'(l) = 0.577215 • • •, and the prime ' denotes differ- 
entiation. This will be proved in Section ^j, equation (35). 



With (p5[), expression (|23| ) can then be analytic continued from the part of 
Qhor-S^grid lying i n the log region into that lying in the lower region. 

Thus, in Qhorz^grid, we have r G Z~, (cr — r) Z, and the analytic 
continuation as 



flf x r = lim r ( 1 + r + £ ) n(«r-r) lo x 
e^o r(e) 

= ! m „Vf(iTH( logI " wl+r " ff) + 7) ) (26) 

3.2.3 Zero region 

As pointed out above, there are two possible choices in this region from the two 
views or schools of thoughts: 

• Type I Fractional Calculus 
Postulate: 

Fractional derivative is abstract. Fractional derivative of a constant can 
be non-zero. So take the Riemann-Liouville Integral as fundamental and 
derive the fractional derivative in this zero region from it. 

Then, 

T(l + r) 

D°x T = — r in the zero region. (27) 

x r(l+r-cr) v ; 

• Type II Fractional Calculus 
Postulate: 

Ordinary derivative of a constant is zero and ordinary derivatives com- 
mute. Fractional derivative should inherit these property from ordinary 
derivative as well — fractional derivative of a constant is zero and frac- 
tional derivates commute. 

D% c = , c is an arbitrary constant 

D° c = D^- 1] [d^ c) = D^- 1] = 0, for c > 1 (28) 

As 

D r x x T = T(l + r) for r > 



10 



by continuity, 

D°x r = for j>r,r>0 
=>• D°x r = in the zero region and commutativity is preserved. 

Hence, in both Types, there is a mutual trade-off. 

In Type I, we chose analyticity and lose commutativity in the zero region. 

In Type II, we chose to preserve and carry over commutativity and lose ana- 
lyticity at the edges of the zero region. 



3.2.4 Entire real (a, r) plane 

The combined analytic continuation is then 
• Type I 

T(l+r + e) x^~^ 



lim 



o r(e) r(l+r-<7) 



log x - (%fj(l+r-a) + 7) 



D£x r = { 



in QhnrA ^ 



horzX^ tgrid 



r (1+r+e) ^ 

I e^o r(l+r+e-cr) 

Type II (fractional derivatives commute) 




elsewhere 



(29) 



in the zero region 



DxX r = < 



r(l+r+e) x^ ( 
& T{e) T(l+r-a) ( logX " ^ l+r ~^ + l) 



T(l+r+e) 
hm — ; - 

e->o T(l+r+e-a) 



x 



in horz\^ grid 



elsewhere 



(30) 



3.3 Analytic continuation from real a to complex s 

To analytic continue D a to D s where complex s = a+it, we just generalise a 



to s in the combined expressions ( |29|) and (29) for for D£x w with complex s, uu. 

For Type II Fractional Calculus, commutativity is preserved. The differ- 
ential operator d/dx commutes with itself and with its inverse, the integral 
operator J dx, and so do D s . 

So in Type II, the commutative operator D s splits into two, each acting on 
functions independently. 

D s _ D (a+it) = D a D it _ D it D a (3^ 
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as illustrated in the commutative diagram of Figure ^, with all the limits, if any, 
taken only at the end after the actions of all the operators have been performed. 



D u 



D" 




Figure 3: Commutative arrows in (s,r) diagram of D£x r in Type II. 

Consider D£x w with complex s = a+it , w = u+iv. 
In the zero region of (a, u) plane, 

D%x u =0 => D^D^x u = D°x u = 
Similarly, in the zero region of {it, iv) plane, 

D it x iv _ q ^ D^D^x iv = D*x iv = 

We can think of the above as D lt expanding the triangular zero region 
of the (a, u) plane into a wedge-shaped volume of infinite length along the t- 
direction, and D a similarly expanding that in (it, iv) plane into another along 
the a-direction as shown in Figure |j. 

This defines the zero space in which D£x w = . 




wedges of zero spac< 




Figure 4: The zero space wedges in (s, w) diagrams where D£x w in Type II. 



For the line intervals in Qhorz\Qgrid > T) lt extends these intervals in the (a, u) 
plane along the t-direction to planar sections, and similarly for D° . 

By now, we have completed the task of defining the action and obtained a 
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complete picture of the operator D s . In essence, what we have done is analytic 
continuing an operator to act on a function complex s-times. Let us carry 
on with applications of the analytic continued D s to find out more about its 
properties and consequences. 



4 Application: Number Theory — 

Generating Integrals and 

Analytic continuation of Finite Series 

"Nature laughs at the difficulties of integration. " — Laplace 

4.1 Generating Integral of Finite Harmonic Series 

It is well known lllf that finite harmonic series 



h(n) = £_ = ! + _ + ... + _ 

= V(l + ™)+7 (32) 
= log n + 7 + 0(l/n) 

Now, we observe that finite harmonic series h(n) appears as the coefficient 
of x n term when we repeatedly integrate logx. 



X 







log x (dx) = x( log x — 1 ) 



/ logx (dx) 2 = ^-(logx- -) 
Jo 11 



/ \ogx(dx) n = — ( log x — h(n) ) 
Jo n\ 

So by analogy to the concept of generating functions, we take 

logx (dx) n as the generating integral of finite harmonic series h(n) 



and so the natural analytic continuation of the generating integral takes the 
form of 

X 

logx (dx) p = ———(logx- h(p)) (33) 
r(l+p) 



From (12^ , 

log x = lim — (x e — 1) 
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The analytic continuation of the integral on l.h.s. of 

as 



can then be evaluated 



PX p X 

/ \ogx(dx) p = lim - / (x € - l)(dx) f 
Jo e Jo 



lim — 

e—0 e 

lim — 

x f 

lim — 



o 

r(l+e) x e 
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Din 



l(dx) 



1 



e^o e |T(l+e+p) r(l+p) 



(34) 



Equating this to the expression on r.h.s. of (|33| ) gives the analytic continu- 
ation of finite harmonic series 



hip) = log a; — lim - 



log x — lim — 



r(i+6)r(i+ P ) 

r(i+e+p) 
' r(i+e)r(i+ P ) 



x 



1 



+ 



r(i+e+p) 
r(i+e)r(i+ P ) 



log a; — lim - (x e — 1) 

<=-+0 e 



r(i+e+p) 
r(i+e)r(i+ P ) 



r(i+e)r(i+ P ) 
r(i+e+ P ) 

- i 



+ lim — 

e^o e 



r(i+e+p) 
r(i+ e )r(i+p) - 

T(l+e+p) . 



1 



logx — lim — (x e — 1) 

e^o e 



+ lim — 

e^o e 



r(i+6)r(i+ P ) 
r(i+e+ P ) 



lim — 



r(i+e)r(i+ P ) 



lim -r(l+p) 



r(i+e+p) 

r'(l+e+p) r(l+e) - r'(l+e) r(l+e+p) 



r(l+e+p) 2 



r'(i+p) 



r(i+p) 



r'(i) 



(35) 



where the limit e — > has been taken with L'Hospital rule. 

([35|) is a result that has been proved in conventional ways in Number Theory. 
It is somewhat surprising and miraculous that it is possible too to rederive it 
with the method of the analytic continued operator D s as above. 



4.2 Series of the Riemann Zeta Function up to Finite Terms 

The analytic continuation (|35|) can be generalised. 
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h(p) 1.5 




2 4 6 8 10 
P 



Figure 5: The curve h{p) = + p)+7 passes through points h{n) = ^ — . 



The Riemann zeta function [12| is defined as 



1 



(36) 



k=l 



The series of the Riemann zeta function up to finite terms is defined as the 
Stegan- Riemann zeta function [11] 



C(s,n) 



n -, oo -. 

E F = coo- E t-s . Re ( fi )> 1 



fe=i 



fc=n+l 

io g (r(i + s)) " 

s=0 



(37) 
(38) 



Its analytic continuation is then given by simply replacing — — with the oper- 

dx n 

ator D w , w G C, 



C(5,«;) = Df log(r(l + s)) 



s=s 
s=0 



(39) 



which can be evaluated when log(r(l+s)) is expressed in the form of asymptotic 
scries 



4.3 Tables of Analytic continued Integrals 

An interesting consequence is that perhaps new editions of Tables of Integrals 
may have to be compiled, eg. compute the coefficient of x k term, w(p,r,a,k), 
in the evaluation of 

' x r {\ogx) a {dx) p 
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5 Application: Calculus — 

Non-integer Power Series, 

Breaking of Leibniz rule and Chain rule 

5.1 Non-integer Power Series 

In Type I Fractional Calculus |0j, we would write the power series for exp(x) 
as 

oo 1 1 
exp(x) = lim V — x n+e = Y — x n (40) 

n=-oo ^ > n( J_ \ I 

In Type II, we would write it simply as 

oo -■ 

^W^EroTw 1 " (41) 

since D£x r = for a > r , r > 0. 

Thus, in Type II, operator D acting on a power series real cr-times 

oo oo r/i,n 

D° y a k x k = y a k . K 7 1 . x k -° for a > (42) 

fro fc ff CTl *r(i+*-*) 

where denotes taking the integer ceiling. 

When a is not an integer, r.h.s. of (|42| ) is a non-integer power series. 
Define the notation 

f(a,x) = DZf(x) (43) 

Think of a in the following way: the one- variable function f(x) is extended to 
a two- variable function f(a,x) in which a has now become a variable of the 
extended function. 



00 (-11 

co S (a,x) = D-cos(x) = D£ £ ^— i 



™2fc 



,=o ( 2 *) ! 

(44) 

£ , r(1+ f\^ 



where Wc(<r) 



(45) 



oo C_iU 

\ j y a w x (2fc+l)! 



f HILJML^I- for ff >0 



(46) 
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where 



Ws{a) 



a 1 

2 + 2 



(47) 



exp(cr, x) = i^"exp(x) = flf /J 



fc=0 



k\ 



~ (-i)* r(i+fc) 

^ fc! r(l+fc-<r) 



(48) 



for cr > 



Table 2: Some tabulated values of Wc(a) and Ws(a) 



a 


0.0 


0.5 


1.0 


1.5 


2.0 


2.5 


3.0 


3.5 


4.0 


4.5 




Wc(a) 





1 


1 


1 


1 


2 


2 


2 


2 


3 




Ws{a) 











1 


1 


1 


1 


2 


2 


2 





We then find that 



cos(<7, 0) 



if a G 2 Z 
otherwise 

exp(<r, 0) = 



sin(<r, 0) 



if a £ 2Z+1 
otherwise 



if a G Z 
otherwise 



cos(cr, x) 
sin(<7, x) 
exp(cr, x) 



± sin(o"±l, x) 
^ cos((7±l, x) 
exp(cr±l, x) 



(49) 



cos(<7±2, x) 
sin(cr±2, x) 



in agreement with the definitions of cos(x), sin(x) and exp(x) when a G Z . 
In addition, it can be observed from Figures [6| and |^ that there exist asymp- 



totic limits 



cos(<7, x) ~ cos(x+ — a) 
sin(<7,x) ~ sin(x + — a) 



> V a as x 



oo 



(50) 



exp(a, x) ~ exp(x) 

which remain to be proved analytically. 

In fact, any function with a power series definition, eg. Bessel functions, 
Fourier series, etc., can similarly have a non-integer power series generalisation. 
In addition, the generalisation of cos(x) to cos(cr, x) can be further extended 
to the case of complex s, eg. cos(s,x) = D£ cos(x). 
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cos(0.01,x) 
i 



cos(0.25,x) 




cos(0.5,x 
1 



0.5 


-0.5 



2ji 471 6n 8ii 2 % 4 7C 6 7C 8 71 2 7t 4 71 6 7C 8 71 2 71 4 71 6 7C 8 7[ 







1 


1 



cos(0.75, x 
i 







t 


1 



cos(1.25,x) 



cos(1.5,x 




2ji 4ti 6ti 8ti 2 7t 4ti 6ti 8 7[ 2 7t 4 71 6 7C 8 7C 



27t 4ti 6 71 8ti 

X 

7T 



Figure 6: Asymptotic limit: cos(<r, x) = D x cos{x) ~ cos(x+ — a) as x — > oo. 



5.2 Breaking of Leibniz rule and Chain rule 

By definition, .D™ with integer n obeys Leibniz rule 

D?{f(x)g(x)} = £ ( fj (D x n - k f(x)) (D x k g(x) 

k=0 \ ) 

and D 1 obeys Chain rule 

Dlg(f(x)) = (D}g) (D l x f) 



(51) 



(52) 



but D s with complex s does not in general. 
However, by observing that 

Dl{fg} = (Dlf)D}{fg} + (D* g) Dg{fg] = {D l J)g + {D l x g)f 

=> Dl = (D x f) Dj + {pi g) D] 

we can express D s where s 6 C in terms of nested sums of Z^'s which we can 
evaluate. 



D* = {{Dlf)D} + {Dlg)Dl) S 

\-(\-{(Dlf)Dj + (Dlg)Dl 



1 + E 



fe=i 



A:! 



! { (s-m+1) 



.m=l 



A 1 /) £>/ 1 + (A 1 5) ^ (53) 
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^exp(<7, x) — exp(x) 




5 10 15 20 

X 



Figure 7: Asymptotic limit: exp(cr, x) — exp(x) — ► as i-> oo. 

Now, we can evaluate expressions of the form of D£ {f(x) g(x)} and g(^f{x) 
by simply by substituting D s with the series on the r.h.s. of (|53|). See Section 
U equation (|7l| ) for the problem of convergence of the series. 

6 Application: Group Theory — 

Analytic continuation of Groups: R^ moc [ n ) Groups, 
Pseudo-Groups, and Symmetry Breaking/Deforming in Groups 

6 - 1 R (mod n) Grou P s 

The differential operator, and its inverse — integral operator, can act on dif- 
ferent functional spaces to generate different discrete groups. These are groups 
of operators, ie. groups with operators as elements. 



—f(x) = f(x) 



Order 


Functional Space 


Symmetry Group 


n = 1 
n = 2 

n = 4 


/(*) = 

m = < 
m = < 


exp(x) 

cosh(x) 
sinh(x) 

± cos(x) 
± sin(x) 


Zi = { 


d , d u 

^ ' dx ' ^ <ix 

' <ix } ' dx 2 

id d ^ ^ \ ^ ~ id 

' <ix ' dx 2 ' (ix 3 J ' <ix 4 



If elements of the functional space are extended from functions /(x) to their 
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A cos(x) d_ 

dx ^ dx 

sin(x) -sin(x) 

X t ^ ^ d 

-cos(x) f- 

dx v ' dx 



Figure 8: Z4 group flow diagram of 

— acting on functional space {± cos(x), ± sin(x)}. 
dx 



analytic continuations f(a,x) = D a f(x) with real a, operators D a acting on 
these extended functional spaces will generate continuous groups or Lie groups, 
eg., 

D£ acting on the functional space jcos(er, x) a G [0,4), x € [0, oo)| gen- 
erates a natural analytic continuation of the Z4 group, 



{^'| CT 'G[0,4)},A 4 ^id 



By analogy to the concept of (mod n) congruence in Number Theory, we denote 
this analytic continued group R( mod 4 ) • 



Dff(a,x) = f(a,x) 



Order 


Functional Space 


Symmetry Group 


<t' = 1 


f{o;x) = 


8xp (a, x) 


R (mod 1) - { D x 


a' e[0,l)},l}j 


= id 


a' = 2 


f(a,x) = < 


cosh(cr, x) 
sinh(cr, x) 


R (mod 2) = { D x 


<x'e[0,2)},^ 2 


= id 


a' = 4 


f(a,x) = < 


± COS ((7, x) 
± sin(<r, x) 


R (mod 4) = { D x 


a' e[0,4)},D£ 


= id 



cos(a, x) 



D, 



4-CT' 




IK 



cos(o'+ o,x) 



Figure 9: R^ moc j 4) group flow diagram of 
D£' acting on functional space jcos(<7, x) a £ [0,4)|. 
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For complex s = a+it , s' = a'+it', 
D s ' acting on functional space |cos(s,x) a £ [0,4) , t £ R, x £ [0, oo)| 
generates a Lie group R^ mCK j ^ X R since D lt ' commutes with D a ' and so D' d ' 
acts independently from D a . 

In general, the topology of such analytic continued groups progresses from 
sets of points on a circle S l for discrete groups generated by d/dx, to a circle 
S 1 for Lie groups generated by D£ , and to a 2-dimensional cylinder S 1 XR for 



Lie group generated by D£ +lt as illustrated in Figure 10 



I — I 
l dx> 



^4 
-• 



Group 



R 



(mod 4) 



(0°i 



R 



(mod 4) 



R (mod4) X R 



{D? +if } 



R (mod4) XR 



Topology 



— ► S — > S X R 



Figure 10: Topology change of groups. 



6.2 Pseudo-Groups in Type II Fractional Calculus 

Consider the analytic continuation of the group elements of SO (2) (the group 
of rotation in a plane) in Type II Fractional Calculus, 



R{6) 



cos v sm < 
- sin 9 cos 



where 9 £ [0, 2 vr) 



R(a, 



cos(<7, 9) sin(cj, 9) 
- sin(cr, 9) cos(cr, 9) 



where 9 £ [0, oo) 

(54) 

R(a, 9) forms a set of sets, paramatrized at 2 levels. The set of sets is parametrized 
by a, and each of these sets is further parametrized by 9. Denote the set of 
these sets as SO(2;a, 9). 

Since 

i?(0, 9) £ SO{2) V 9 and R(a, 9) ~ R(0, 9+^ a) as 9 -» oo (55) 

we are motivated to introduce the concept of pseudo-groups. 

A pseudo-group G(p\, p2, ■ ■ ■ , Pk) of a group G is a set which gradually 
acquires the group properties or satisfies the group axioms of G as some of the 
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parameters p\ , p2 , • • • , Pk of the set approach limiting values or tend asymptot- 
ically to infinity. 

50(2; a, 9) is a pseudo-group of 50(2) since it is isomorphic to 50(2) when 

• the parameter tr-nieZ while 9 varies freely in the interval [0, oo), 

R(a, 9) . R(a', 9) ~ R{a+o', 9) as a, u'^n.n'eZ 

lim SO(2;a,9) 50(2) (56) 

• (alternatively) the parameter 9 — > oo while cr varies freely in the interval 
(0,2), 

R(a,9).R{a,9') ~ i?(0, 0+|<r) . i?(0, 0'+|<t) = fl(0, 0+0'+7r<r) 



as 9,9' ^ oo 
lim 50(2; <r, 0) 9* 



50(2) 



(57) 



o i 




S0(2) 



Figure 11: SO(2;a, 9) plane diagram. 



We can define a group property deviation measure W(G(a,x),G \ a,x) for 
a pseudo-group G(a,x), a measure of how much group property the pseudo- 
group has lost or deviated from the associated "parent" group G from which 
it is analytic continued. When the pseudo-group becomes isomorphic to the 
parent group for certain values of the parameter, the measure should be zero. 

For the case of pseudo-group 50(2; a, 9), 
W(SO(2;a,9),SO(2) \ a, 9) 

7T 



R(a,9) -R(9+-a) 



7T 



cos(<7, 9) — cos(0 + — cr) 



7T 



sin(er, 9) — sin(0 + — a] 



sin(<7, 
cosfo", 



erin(0+~<7) 

in 71 \ 

cos{9+—a) 
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Similarly for the simple case of U(l;a,ix), a pseudo-group of £7(1) where 
x e R, 

exp(cr,zxi) exp(cr / ,ix2) ~ exp(cr, i{x\+X2)) as a, a' — ► n, n' G Z 

lim £7(1; a, x) 9* 17(1) 

(T— >nGZ 

exp(<T,ixi) exp(o",ix2) ~ exp(<r, +£2)) as xi,x%— >• 00 
lim Z7(l;er,a?) = £7(1) 

W(£7(l; a, zx), £7(1) | a,ix) = exp(a,ix) — exp(zx) 
This measure was plotted in Figure [?]. 
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6.3 S0(2;a, 6) Rotations and Deformations in Type II 



Figure 13 shows the effect of planar rotations and deformations of SO(2;a,6) 
on a square with vertices { (1,-1), (1, 1), (-1, 1), (-1,-1)} on a sequence of (x,y) 
planes clipped by square windows of size x£ [—2,2], y£ [—2,2]. The deforma- 
tion effects seem to be a combination of rotations and contractions/dilations. 



a\e 



tt/16 tt/8 tt/4 tt/2 3tt/4 tt 3tt/2 2tt 4tt 6tt 8tt 



Figure 13: SO(2;a, 9) rotations and deformations (color illustration). 



6.4 Symmetry Breaking/Deforming in Groups in Type II 



In the Higgs mechanism of Spontaneous Symmetry Breaking [13], 



the symmetry of the effective potential V e ff in a Lagrangian density C 
with respect to a gauge group G is preserved, while 
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• the symmetry of the quantum state ip satisfying the equations of motion 
derived from C is broken and reduced to that of a subgroup, H C G . 

The profile of V e ff changes with energy or temperature. At high energy or 
temperature, the symmetry of ip is restored from H — > G . 

In the case here, the symmetry breaking is very different. The symmetry 
in a group G itself is broken to a subgroup H C G or "deformed" into an 
approximate symmetry of G . 

Take 50(3; (o"i, x{), (<7 2 , #2)), a pseudo-group of SO(3), as an example. 
When both o\,G2 = 0, SO(3; (a±, x±), (cr 2 , £2)) is isomorphic to .SO^). Now, 
choose o"i Z and 02 = . 

{broken | f small ) 

approximate > for < intermediate > x, and the 
restored I 1 large I 

S"0(3) symmetry in a sphere is "deformed" to an approximate SO (3) symmetry 
or completely broken to SO (2) in a plane depending on the chosen values of 
o\ and x. Now, set both a\,a2 Z. The SO(2) is further "deformed" to an 
approximate SO (2) symmetry or broken to Identity. 

The symmetry breaking/deforming sequence is then 

0-1,(72=0 CTl 6? Z, (72=0 (7l,CT2^Z 

50(3;(<7i,a;i),(<T2,X2)) -» ► 50(3) -» > 50(2) -> ► 1 

Similarly for SU(N; (cr l5 xi), (o"2, X2), • • • , (cjv, £iv))> the symmetry break- 
ing/deforming sequence is 

SU(N; (ai,x{), (<t 2 ,x 2 ), • • • , (cr N ,x N )) 

(7l,(72>---,O- J v=0 

^ > SU(N) 



(J\4. Z , CT2,---, (7JV= 

> SU(N-1) 



°"l,---,CTjV-20Z,CTjV_l,<7jV=O 

> 5C/(2) 

C1,---,O-JV-10Z,(7JV=O 

^ > C/(l) 

(71, <T2,-", crjv^ Z 
► 1 

Perhaps this mode of symmetry breaking/deforming in groups might have 
some useful applications for models in Particle Physics and Cosmology. 

In Particle Physics, the symmetry of the flavor of quarks are not exact 
symmetry but only approximate symmetry of Gell-Mann's Eightfold way SU (3) 



1 14] or GUT SU (5) []15[ because different flavors of quarks have different masses. 
Light quarks do not transform exactly into heavy quarks, and perhaps the 
effects from the presence of gluons and glueballs in composite particles need to 
be added into the symmetry. Perhaps SU(3) and SU(5) can be "deformed" in 
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this way to an approximate symmetry that will fit the phenomenological data 
better. 

Now, take a pseudo-group to describe the product of residual exact and ap- 
proximate symmetries of present day Universe. If we set x oc T, the temperature 
of the Universe, as we go back in time, the temperature T goes up, x goes up, 
and we find that the approximate and other broken symmetries are gradually 
being restored. The rate the symmetries are being restored will be dependent 
on the values of a%, 02, • • • , o-jy, the parameters of the pseudo-group. The fully 
restored symmetry will be the symmetry of the parent group of the pseudo- 
group. Qualitatively, this model resembles the unification of gauge groups in 
Cosmology |l6|]. It might be interesting to study and develop this mode of 
symmetry breaking/deforming for approximate symmetry groups (eg. iso-spin 
group, Eightfold way SU(3), and GUT SU(5) ) as well as for the gauge groups. 



7 Application: Algebra — 

Analytic continuation of Dirac Equation and Algebra 

"I think that there is a moral to this story, namely that it is more 
important to have beauty in one's equations than to have them fit 
experiment." — P.A.M. Dirac in Scientific American, May (1963) 

Dirac equation Klein-Gordon equation 

i —ip = (^— ia . V + f3 rnj ip =>■ (i ■—) ip = (^—ia . V + (5 rnj ip 



&t T \ ^ r V dtj 

dt 2 

ip = (— V 2 + m 2 *\ip where ^ = ^ ^ • • • ^ 



rip = -^ctiCtjViVj - im(^af3 + f3aj.V + f3 2 m 2 \ip 



'■j 



dt 2 

giving Dirac Algebra |1 



{a i ,a j } = 2 5 ij l , {a* , 0) = , af = (3 2 = 1 (59) 
f d 1 

Think of < i — , — i V > as a basis. The basis can be analytic continued with 



dt 

the D s operator. 

Dirac equation can then be analytic continued to 

e n i/p D 2/p^ = i/p g(2 /p) _ R (2/p) + p{2/p) m 2/p^ ^ (6Q) 

where p = 1 , fl60| ) Klein-Gordon equation 
p = 2 , ( pC|) =>■ Dirac equation 
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Introduce the notation for generalised symmetrisation 



^perrr^ ^permj • • • ^perm p (61) 



perm(ii,i 2 ,...,j p ) 

where the sum is over all permutations of the p indices. 

The generalised symmetrisation can be re-expressed in terms of a sum of 
permutations of nested anti-commutators, eg. 

{«, b, c} = ({{a, b},c} + {{b, c}, a} + {{c, a}, b}) (62) 

{a, b, c, d} = 

1 / {{{a,b},c},d}+{{{b,c},d},a}+{{{c,d},a},b}+{{{d,a},b},c}\ 
+ {{{a, 6}, d}, c} + {{{6, c}, a}, d} + {{{c, d}, b}, a} + {{{d, a}, c}, 6} 
\+ {{{a, c}, 6}, d} + {{{a, c}, d}, b} + {{{6, d}, c}, a} + {{{6, d}, a}, c} / 

(63) 

Now, for p = 3, 
L.h.s. : 



2(2!) 



R.h.s. : 

(e~™/3 aW.D 2/z + /? (2/3) m 2 / 3 )% 



52_ 

dt 2 



i,j,k 



(2/3)^(2/3)^(2/3)^2/3^2/3^2/3 



+ e- 27ri / 3 m 2 / 3 ^ 



/ af /3) af 3) /3 \ 

+ a? ^ af s > 
V+/3 af /3) af /3) 



D 2/3 D 2/3 



+ e--/ 3 m 4 / 3 



a 



(2/3) 



£(2/3) 



+ £(2/3)^,(2/3) £(2/3) 



+ 



£(2/3) 



(2/3) 



J 



+ 



m 



[£(2/3)]' 

+ e- 2 -/ 3 m 2 / 3 i^{af /3) D 2/3 , af 3) D 2/3 , /?} 



+ 



+ 



e --/3 m 4/3 |^(2/3) } £ )/3 j_ £(2/3) 



£(2/3)" 



m 
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o 



(2/3) 



.D 2 



{af/ 3 ),af 3 ),4 2 / 3 )} = 3!^l , 

{af /3) , af /3) , /3( 2 / 3 ) } = , {af /3) , /?( 2 / 3 ) , /3( 2 / 3 ) } = 



a 



(2/3) 



0(2/3) 



For general p, Dirac Algebra is analytic continued to 

ia {2/p) a (2/p) a {2/p) a (2/p) a (2/p) \ - o' <J ■ ■ ■ 1 



(ag/rt, ^),. ^)) = , 
(p-i)-times 



(2/P) 



^(2/p) 



(64) 



From here, we may proceed on to find representations of this analytic con- 
tinued algebra and study the properties of the associated analytic continued 
spinors. Perhaps they have interesting properties. 

As a hint, even the matrix representation of the finite difference of D s itself 
has surprising properties. To this we turn to next. 

8 Application: Matrix Representation — 

Analytic continuation of Matrices, 

from local Finite Difference to non-local Finite Difference 



"We [he and HalmosJ share a philosophy about linear algebra: we 
think basis-free, we write basis-free, but when the chips are down we 
close the office door and compute with matrices like fury. " — Irving 
Kaplansky in Paul Halmos: Celebrating 50 Years of Mathematics 

Given a generic matrix [M], we all know how to compute [M] n , the matrix 
[M] raised to an integer power n G Z. It is just trivially multiplying the matrix 
[M] by itself n-times. 

Now we wish to compute [M] CT , the matrix [M\ raised to a real non-integer 
power a G R. A generic matrix may have degenerate eigenvalues and so cannot 
in general be diagonalized. However, we can obtain [M]°~ as follows: 

For rational a = p/q, where p, q G Z, 
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[TVf ]p/<? ca n be obtained by solving for each element of the matrix [A] in 
the matrix equation [A} q = [M] p 

since formally ([M] p / q ) q = [Mf' q [M] p / q ■ ■ ■ [M} p ' q = [M] p = [A] q 



and so [A] = [Mf/ q . 



g-times 



In Finite Difference, if we choose the matrix representations of differentia- 
tion D 1 to be 

/ 1 -1 \ 
1 -1 



[ D x] = 



1 -1 

1 / 



'(Ax) 



(65) 



with 



/ 1 -2 1 

1 -2 1 



m = [D l x ] [££] = [R 



v 



1 -2 1 
1 -2 

1 / 



(Ax) 2 (66) 



and that of integration D 1 as the inverse of D 1 , 

/ 1 1 1 1 



ID' 1 } = [A 



X J 



1 1 1 



1 1 

1 / 



(Ax) 



(67) 



we then have 



m-timcs 



\m = \dt = n\ • • • \Di 
[D- m ] = [D i x \- m = \di\-* m-i-m 



X J 



/ 1 -1 

1 -1 



[ D x\ [/(*)] 



rrt-times 
\ ( fir,,) \ 
f{x n -i) 



V 



1 -1 
1 J 



f(x 2 ) 



(Ax) 



/ 1 1 1 1 
1 1 1 



ID' 1 ] [f(x)} = 



1 1 



\( f(Xn) \ 




f(x n -i) 






(Ax) = 


f(X2) 




A f(xi) J 





( Dxf(Xn) \ 

D*f{x 2 ) 
\f( Xl )/(Ax)J 

D-'fix^f 
D-'fix)^- 1 



D^f{x)\ x x \ 
\ f(x n )(Ax)J 
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where Xk = x\ + (k — 1) Ax. 

[D^] have degenerate eigenvalues and thus cannot be diagonalized. How- 
ever, following the above approach, we can compute the matrix representation 
of D" a E R. It can be verified that 



( K<*> 1) 


b(a,2) 


b(a, 


3) •• 


b(a,n) \ 












b(a,l) 


b(a, 


2) •• 


6((7, n-1) 














b(a, 


1) •• 


6((7, n-2) 


1 (Ax) a 


for 


(7 > 


V 








6((7,1) ) 










/ b(-a, 1) 


b(-a, 2) 


K 


a,3) 


■ ■ ■ b(-a, n) 


\ 


-1 








b(-a,l) 


K 


(7,2) 


■ ■ ■ b(-a, n-1 


) 












K- 


(7,1) 


■■■ 6(-<7,n-2 


) 


{Ax)- a 


for 


(7 < 


V 








b(-a,l) 


) 









satisfy 

where b(a, k) 



(_!)(*-!) 



for 



(7 



<7l + 0"2 H h CTj, 



(68) 
(69) 



lim 



r(l+(7+e) 



which is incidentally the /c-th of 



r(ife) e-*or(2+(7+e-A;) 

the binomial expansion (1 + (—1))°". 

The matrix representation [D m ] is sparse while [D a ] is in general dense - 
all the elements in the upper tri-diagonal block become non-zero. 

In Finite Difference, sparse matrix entails taking the differences between 
only neighboring sets of points, while dense matrix entails taking the differences 
among points almost everywhere in the domain — a non-local effect. 

From 



hmjZT] [f(x)] 



lim 5>(a,fc)/( 



%m—k-\-l , 



(Ax) 



(70) 



k=l 



where 1 <C m < n. 

If /(x) is an integer power series, D a f{x m ) on l.h.s. is in general a non- 
integer power series. The corresponding matrix representation on r.h.s. is a 
sum of ordinary integer power series. The non-local effect can be then seen 
to arise from approximating the non-integer power series by a sum of ordinary 
integer power series. 

This is in parallel with the application of fractional derivative as a pseudo- 



differential operator in non-local field theory by Barci et al [18]. 

For the case of [M] u and [M] s , the matrix [M] raised to irrational u and 
complex s respectively, we turn to the series expansion method in equation (fn|) 
below. 
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9 Analytic continuation of Generic Operators 



Now, let's go beyond the analytic continuation of differential and integral op- 
erators to the analytic continuation of generic operators. 

We are used to thinking of an operator acting once, twice, three times, and 
so on. However, an operator acting integer times can be analytic continued to 
an operator acting complex times by making the following observation: 

A generic operator A acting complex s-times can be formally expanded into 
a series as 
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(71) 



n-l 

ik- 

k=0 

n-l 

is the identity operator. 

In the nested series on r.h.s, all the operators A's are raised to integer powers 
which we can evaluate as usual. The region of convergence in s and the rate of 
convergence of the series will in general be dependent on operator A, parameter 
w, and the operand on which A acts. The resulting series then defines A s , the 
analytic continuation of the operator A, in the region of s where it converges. 



10 Problems and Challenges 

"Mathematics is not yet ready for such problems," — Paul Erdos 
in The American Mathematical Monthly, Nov. (1992) 

10.1 Analytic continuation of Bernoulli numbers and polyno- 
mials, a new formula for the Riemann zeta function, and 
the Phenonmenon of Scattering of zeros 



Examples of interesting mathematical applications are analytic continuation of 
Bernoulli numbers and polynomials, the derivation of a new formula for the 
Riemann zeta function in terms of a nested series of Bernoulli numbers, and 
the observation of particle-physics-like scattering phenomenon in the zeros of 
the analytic continued polynomials as described in S. 



For instance, an operator was found in [19] to generate Bernoulli numbers. 
Applying the series expansion to the operator analytic continues the Bernoulli 
number to a function 
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(72) 
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which was verified to converge, for Re(s) > (1/w) , real w > , to 

B(s) = - s C(l - s) (73) 




Note that B(n) = B n for n > 2 but 5(1) = 1/2 while B 1 = - 1/2. 
It was then realised that there is actually an arbitrariness in the sign convention 
of B\. The analytic continuation of Bernoulli numbers fixes the arbitrary sign 
convention, and requires the generating function of Bernoulli numbers to be 
redefined for consistency as 
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(74) 



(75) 



which only changes the sign in the conventional definition [|ll]] of the only non- 
zero odd Bernoulli number, B%, from B% = — 1/2 to !?i = -B(l) = 1/2 . 



(0), 



From ( 72 ) and (O) , by the functional equation of the Riemann zeta function 
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a nested sum of the Riemann zeta function itself evaluated at negative integers, 
which converges for Re(s) > (1/w) , real w > , with redefined B\ = 1/2, 
and the limit only needs to be taken when s = 1,3, 5, . . . £ the set of 

positive odd integers, for which the denominator cos I — I =0. 



10.2 Quantum Operators 

"The modern physicist is a quantum theorist on Monday, Wednes- 
day, and Friday and a student of gravitational relativity theory on 
Tuesday, Thursday, and Saturday. On Sunday he is neither, but is 
praying to . . . find the reconciliation between the two views. " 
— Norbert Wiener 



Quantum Mechanics, Quantum Field Theories and Canonical Quantum 
Gravity are full of non-commutative operators. 

As a start, let's try to apply the idea of analytic continuation of operators to 
the creation and annihilation operators of a simple harmonic quantum oscillator 

m 

a) \n) = \/n+l |n+l), a\n) - 
From the nested series expansion (|7l|), 



n \n — 1) , a |0) 







(77) 



(wl — [wl — a] 




m=l 



W 



P 

,m. 



and similarly for cd^y 

Consider a s acting on the quantum state \n), 
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(78) 
(79) 



The nested series expression is straightforwardly computable but the inter- 
pretation is not clear. 

When s is not an integer, the inner series nested within the outer series on 
r.h.s. of (ff^) terminates at the (p = n)-th term. However, the outer series does 
not terminate but is an infinite sum for both operators. On the other hand, 
the eigenvalues of the creation operator for simple harmonic quantum oscillator 
are not bounded above, under the action of a\ S s, In) seems to transform into a 

(w) ' I ' 

superposition of \n) and an infinite towers of all other states higher than |n). 



33 



Case s = /c G Z + , w £ R or C : 

w m {k,n,w) = 



y/n\/(n — k)\ {m = n—k) 
otherwise 



(80) 



and so we recover (77) from the nested series expansion ([78]) independent of w. 



Case seR, s Z + , w G R or C : 

There may be divergences in the nested series expansions. Much remains to be 
worked out and clarified. 

Now, we turn to the commutators of these analytic continued operators. 

f s f 

A formal nested series expansion of the commutators of and a\ m ^ is 
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(81) 



From the canonical commutation relations pjj , 

[a, a f ] = 1, [a, a] = = [a f , a f ] (82) 

the [a m , a) 17 *] in the nested series can be evaluated as usual, and thus the nested 
series expansion (Blj) is formally computable. 

Similar generalization applies to fermionic operators satisfying Grassmann 
algebra, SUSY operators 



{Q,Qt } = 2H/u^{Ql w) ,Q^ l) }, 



Virasoro generators |22j in String theories 



(n-mjLn^m + — n(n 2 -l) 5. 
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(83) 



(84) 



its jd 
[■ u n(w) i ±J m{w' 



Superconformal algebra [23] in Superstring, their respective vertex operators, 
and Lie algebra in general. 
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These generalizations seem to have interesting mathematical structures. 
Further aspects and detailed computations will be presented and explored in a 
forthcoming paper [24]. 



11 Conclusion 

Analytic continued operators have been demonstrated to exhibit intriguing 
properties. In addition, fractional derivatives in the conventional Riemann- 
Liouville Fractional Calculus do not generally commute but an extension in 
which they commute has been found and applied to various fields. These meth- 
ods of analytic continuation of operators may after all turn out to be a general 
and powerful exploration tools in Maths, Physics, Sciences, and Engineering. 

Calculus is never quite the same again. It would be interesting to imag- 
ine what Newton and Leibniz would say on this analytic continuation of their 
discoveries — Calculus, and Dirac of his equation - Dirac equation. 

Perhaps the most unexpected, and yet "inconsequential", consequence is 
that Figure [14] clearly points out that the commonly adopted definition of the 1st 
Bernoulli number B\ has the wrong sign. There was actually arbitrariness in 
its sign convention and the analytic continuation of the operator that generates 
Bernoulli numbers || fixes that arbitrariness, requiring that B\ = — 1/2 to 
be redefined as B\ = 1/2 for consistency. However, the B\ = 1/2 definition 
has been so widely used — in every Math, Physics and Engineering book or 
paper where Bernoulli numbers appear, one almost certainly find B\ = —1/2. 
I can only hope that the readers will be persuaded in the light of this new 
mathematical fact to change and adopt the consistent definition B\ = 1/2 and 
the corresponding defining equations ([74|) and (|75[). 
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